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A higher order QED calculation of the ultraperipheral heavy ion cross section for \x~ pair 
production at RHIC and LHC is carried out. The so-called "Coulomb corrections" lead to an even 
greater percentage decrease of fj> + fj,~ production from perturbation theory than the corresponding 
decrease for e + e~ pair production. Unlike the e + e~ case, the finite charge distribution of the ions 
(form factor) and the necessary subtraction of impact parameters with matter overlap are significant 
effects in calculation an observable ultraperipheral /i + fx~ total cross section. 
PACS: 25.75.-q, 34.90. +q 



I. INTRODUCTION 



In recent years sufficient progress has been made in 
evaluating higher order QED corrections to e + e~ pah- 
production in ultraperipheral heavy ion collisions to al- 
low a meaningful comparison with data. The comparison 
of calculated e + e~ pair production rates 1] with STAR 
data0 provides the first evidence of higher order QED 
effects at RHIC. The analogous higher order corrections 
to pair production are now of interest with the 

anticipated program of ultraperipheral heavy ion colli- 
sions at the LHC[3|. In particular it has been suggested 
by Kharlov and Sadovsky (see [3|) that since muon pairs 
are easy to detect (in the ALICE detector) and simple 
to calculate (in perturbation theory) they can be used 
as a luminosity monitor at the LHC. And a recent paper 
by Hencken, Kuraev, and Serbo has presented approx- 
imate analytical calculations indicating that the higher 
order Coulomb corrections to ji + pT pair production are 
smallQ. However, in this paper I will present less ap- 
proximate numerical calculations showing that higher or- 
der QED corrections to the rates are not small (and 
thus would require inclusion e.g. in calculations to obtain 
a meaningful luminosity measurement). I will further ar- 
gue that an additional reduction in /x + fj,~ rates, relatively 
insignificant in the analogous e + e~ process, arises from 
unitarity corrections necessitated by the exclusion of ion 
trajectories with nuclear overlap. 



In section II I discuss the method utilized here to calcu- 
late higher order QED lepton pair production with heavy 
ions. Section III considers the differences between rates 
for e + e~ and pair production in the usual impact 

parameter independent calculation: how they scale and 
the relatively larger effect of the heavy ion form factor 
for the pairs. Section IV introduces the impact 

parameter dependent representation necessary for con- 
sideration of unitarity effects, especially the exclusion of 
ion overlap trajectories. A summary of the results for 
computed total /i + /.t~ pair production cross sections at 
RHIC and LHC is presented in Section V. 



II. HIGHER ORDER QED: METHOD OF 
CALCULATION 



The method employed in this paper for higher order 
pair production is an extension of previous work 
on higher order e + e~ pairs. A broad review of issues in- 
volved in the e + e~ pair calculations has been presented 
in a recent review^, and a discussion of the factoriza- 
tion of the different electromagnetic processes and the 
applicability of the semiclassical description is found in 
Ref. [6| . Here I review the history relevant to the present 
calculations. 

The possibility of accurate higher order e + e~ pair cal- 
culation originated with the realization that in an appro- 
priate gauge 0, the electromagnetic potential of a rel- 
ativistic heavy ion is to a very good approximation a 
delta function in the direction of motion of the heavy ion 
times the two dimensional solution of Maxwell's equa- 
tions in the transverse direction 0, [9[. This led to the 
closed form solution of the time-dependent Dirac equa- 
tion for lepton pair productionfiol. [ill [l2T| . However this 
original solution needed to be corrected for the follow- 
ing reason. One apparent consequence of the original 
solution was that rates for pair production in the exact 
solution agreed with the corresponding perturbation the- 
ory result. It was subsequently pointed out by Ivanov, 
Schiller, and Serbo [l3| that this heavy ion conclusion was 
contrary to the well known fact that photoproduction of 
e + e~ pairs on a heavy target shows a negative (Coulomb) 
correction proportional to Z 2 that is well described by 
the Bethe-Maximon theory [3]. These authors went on 
to compute large Coulomb corrections to the pair total 
cross section by considering higher order Feynman dia- 
grams in a leading logarithm approximation. Lee and 
Milstein[l5|, [l6| came to essentially the same result for 
the Coulomb corrections. They pointed out that the orig- 
inal Dirac equation solution involved an integral over the 
transverse spatial coordinates that was not well regular- 
ized. Lee and Milstein constructed an appropriate regu- 
larized transverse integral in the low transverse momen- 
tum (k) approximation that could be solved analytically 
to obtain the Coulomb corrections. They also noted that 
replacing the original transverse potential — 2iZaln(p) 
with 2iZaKo(pu>/j) gives a properly regularized expres- 
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sion for the original transverse integral 



a simple form factor 



F(k) = 2n J dppJo(kp){exp[2iZaK (pij/-f)} - 1}, (1) 
that goes over into the correct lowest order expression 



F (k) = 



k 2 + U! 2 /j 2 



(2) 



in the perturbative limit. The modified Bessel func- 
tion Ko(pu)/j) = — ln(p) plus constants for small p and 
cuts off exponentially at p ~ j/lu, where 7 is the rela- 
tivistic boost of the ion producing the photon and u) is 
the energy of the photon. I previously carried out nu- 
merical calculations utilizing these expressions and ob- 
tained results identical to those of Lee and Milstein in 
their small k limit [171] . In my previous cross section 
calculations [Lsl [l9| and in what follows, Eq. (1) is uti- 
lized for the higher order calculations and Eq. (2) for 
lowest order. 

For impact parameter dependent cross sections the cal- 
culations presented here make use of the methods of cal- 
culating e + e~ pair probabilities previously described [l^]. 
The impact parameter (b) dependent amplitude presents 
a particular numerical challenge since it involves a rapidly 
oscillating phase exp(ik • b) in the integral over the trans- 
verse momentum k transfered from the ion to the lep- 
ton pair. The usual method of evaluating the perturba- 
tive impact parameter dependent probability is to first 
square the amplitude and then integrate over the sum 
and difference of k and k'. Here I have integrated be- 
fore squaring, and I deal with the rapid oscillations with 
the piecewise analytical integration method previously 
described [l^]. In that previous ^-dependent calculation 
of the total cross section for e + e~ production, half of 
the contribution comes from b > 5000 fm and contribu- 
tions up to b = 10 6 fm are considered. Due to the large 
values of b contributing, that calculation was somewhat 
crude. However integration over b reproduced the known 
cross sections calculated with the b independent method 
or calculated from the very accurate analytical Racah 
formula[2(| to about 3%. It can also be noted that the 
computed perturbative b dependent probabilities in that 
paper were in relatively good agreement with calculations 
in the literature 2l| available for b < 7000 fm. 



III. 



SCALING OF p+p- WITH 
SECTIONS 



CROSS 



Let us begin by reviewing the scaling of p ■ p~ cross 
sections from the corresponding e + e~ cross sections. For 
point charge heavy ions (no form factor) if length is ex- 
pressed in terms of 1 /mi and energy in terms of mi then 
the total lepton pair cross section a(p + p~) is identical 
to er(e + e - ) 

A form factor g(k) may be defined that modifies the 
expressions for F{k) in Eqns. (1) and (2). If one assumes 



l + k 2 /A 2 



where for Au or Pb 



A ~ 80 MeV = 160 m e 



.75 m 



then Eq. (2) for the perturbative limit becomes 



{k 2 +L0 2 h 2 ){\ +k 2 /K 2 Y 



(3) 



(4) 



(5) 



k is cut off at the low end when k 2 -C (w/7) 2 . In the 
perturbative case for e + e~ pairs it has been shown that 
the effect of the form factor seems to be present only 
where the impact parameter is of the same size as the 
nuclear radius [12] ■ However the situation is different for 
p pairs. At the high end the form factor cuts off when 
k 2 ^S> A 2 . The form factor contributes if A 2 is comparable 
to or less than (w/7) 2 , the cutoff of k without the form 
factor. Assume that at this high end cutoff without the 
form factor k ~ lOOw/7 ~ 10 for RHIC. Clearly for p 
pair production the sum of the us for the two virtual 
photons must be greater than twice the mass of the muon. 
Thus for even the lowest energy p pairs (corresponding 
to large impact parameters) at least one u) > m^ and the 
form factor is important. On the other hand, the form 
factor is relatively insignificant for the total cr(e + e _ ) and 
contributes only at electron energies some two orders of 
magnitude above the electron mass, comparable to the 
value of A. Without a form factor 



^=(^) 2 = 2,34X10- 5 
cr(e+e ) m^' 



(6) 



But with a form factor the perturbation theory result 
calculations give 

^^l=0.61xl0- 5 = .26x(^) 2 (7) 
cr(e + e ) m^ 

for RHIC, and 

^J^l = 1.16 x HP 5 = .50 x (^L) 2 (8) 
a(e + e ) m^ 

for LHC. 

To include a form factor in the eikonalized expression 
with Coulomb corrections Eq. (1) then the most obvious 
prescription is to apply the form factor to the transverse 
potential: 

F f (k) = 2n J dppJ (kp){exp[2iZag(k)K (pLj/"/)} - 1}. 

(9) 

This expression obviously goes into the correct perturba- 
tive limit Eq. (5). A simpler expression is to take the 
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form factor only to first order but the Coulomb correc- 
tions to higher order 

F f0 (k) = 2irg(k) J dppJ Q (kp){exp[2iZaK (pu/j)] - 1}. 

(10) 

Again, this expression obviously goes to the correct per- 
turbative limit Eq. (5). This is the expression that will 
be used in this paper. A discussion of the validity of this 
approximation is given in Appendix A. 

For simplicity in calculation and simplicity in compar- 
ing with Ref.Jj], the form factor g{k) Eq. (3) has also 
neglected any dependence on longitudinal momentum. 
Including a longitudinal momentum dependence would 
make a small reduction in cross section values, about 5% 
for RHIC and 1% for LHC, as discussed in Appendix B. 

I have previously calculated [3| that there is a 17% 
reduction at RHIC and a 11% reduction at LHC in the 
exact total <r(e + e~) from the perturbation theory result. 
For the er(/i + /i~) here the corresponding reduction from 
perturbation theory is even greater, 22% at RHIC and 
14% LHC. The present perturbative a(p + p~) calcula- 
tions are in fairly good agreement with the calculations 
of Hencken, Kuraev, and Serbo [H, but the present ex- 
act cross section calculations are in disagreement with 
their argument that Coulomb corrections are relatively 
insignificant for p pairs. 

So far the calculations presented have been performed 
in the impact parameter independent representation. 
There is an additional reduction that comes into play for 
an observable a(p + p~) that arises from unitarity consid- 
erations, and one must make use of the impact parameter 
representation discussed in the following section. 



IV. IMPACT PARAMETER AND UNITARITY 

The perturbative (Born) cross section and correspond- 
ing cross sections with higher order Coulomb corrections 
discussed in the previous section correspond to an inclu- 
sive cross section, constructed from a probability corre- 
sponding to the number operator for a given process. If 
one considers an exclusive cross section, e.g. exciting a 
p pair and nothing else in a heavy ion reaction, then one 
must consider unitarity corrections for competing pro- 
cesses in an impact parameter representation as will be 
seen below. For p pair production the main unitarity 
corrections arise in principle from competing e + e~ pair 
production, Coulomb dissociation of the heavy ions, and 
nuclear processes at ion-ion overlap. 

Hencken, Kuraev and Serbo have observed that while 
unitarity corrections are small for e + e~ cross sections 
they are large for corresponding p + p~ pair production 
The perturbative Born cross section for e + e~ produc- 
tion, corresponding to an inclusive cross section, is lit- 
tle increased from the exclusive cross section. The per- 
turbative Born cross section for p + /i~ production also 
corresponds to the inclusive cross section, but the ex- 
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FIG. 1: Unitarity reduction factor for relevant processes in the 
calculation of exclusive \x pair production at RHIC: dashed 
line - nuclear collisions; dotted line - Coulomb excitation of 
one of the heavy ions; solid line - nuclear collision or Coulomb 
excitation; diamonds - e + e~ pair production. 



elusive cross section is significantly reduced by unitar- 
ity corrections due to the simultaneous production of 
e + e~ pairs along with the p + p~ pairs. Lee and Mil- 
stein have recently developed a quasi-analytical proce- 
dure to include the higher order Coulomb corrections in 
calculating the imp act parameter dependence the e + e~ 
pair production[23| . Based on their procedure Jenschura, 
Hencken and Serbo have updated the consideration of the 



e pair unitarity corrections 24 1 . 



In practice some unitarity corrections are relevant to 
what is actually measured and some are not. While it 
is an enlightening theoretical exercise to consider e + e _ 
pair unitarity corrections to p pair rates, in practice the 
dominant contributions of soft e + e _ pairs are of an en- 
ergy scale orders of magnitude too small to be observed 
in an experiment designed to observe p pairs. On the 
other hand, when calculating p pair rates in an impact 
parameter representation, a correction must be made to 
exclude the lowest impact parameters of ion-ion overlap, 
where the dominant processes are nuclear. 

If one assumes independence of the various heavy ion 
reaction processes, then the probability of a reaction 
leading to a final state differing from the incoming chan- 
nel is given by the usual Poisson distribution. If P(b) 
is the sum of the probabilities for producing an excited 
state by a heavy ion reaction at a given impact parameter 
b 



P(b) = HiPi(b), 



(11) 



where each Pi(b) is the inclusive probability of a specific 
final state i, then the exclusive probability for a given 
final state P"(6) is 



i?(6) = Pi(b)exp(-P(b)), 



(12) 
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FIG. 2: Impact parameter dependence of contribution to total 
cross section for p pair production at RHIC: dashed line - 
perturbation theory; solid line - higher order calculation. 



FIG. 3: Impact parameter dependence of contribution to total 
cross section for e + e~ pair production at RHIC: dashed line 
- perturbation theory; solid line - higher order calculation. 



and the probability of remaining in the initial state is 

PZ(b) = l-exp(-P(b)). (13) 

Figure 1 shows the unitarity reduction factor 
exp(~P(b)) in Eq. (12) evaluated for the probabilities 
of various processes for the case of Au + Au at RHIC. 
In agreement with the previously discussed above work 
in the literature [1, [23|, [24[ the unitarity effect of e + e~ 
pair production (diamonds) is significant for low and in- 
termediate impact parameters. Following the methods of 
Ref. [Ill, [2(| I have also calculated the unitarity reduction 
factor exp(—P(b)) for Coulomb dissociation and nuclear 
dissociation. 

It is instructive to compare the impact parameter de- 
pendence of the contribution to fi pair production and 
e + e~ pair production at RHIC. Figure 2 shows the dis- 
tribution for for /i pairs and Figure 3 for e + e~ pairs. 
In both figures dashes correspond to perturbation the- 
ory and the solid line the higher order calculation. The 
shift in scale mentioned in the previous section is evi- 
dent. Comparing the region of e + e~ reduction shown in 
Fig.l with the regions of dominant cross section contri- 
bution for /i pairs (Fig. 2) and e + e _ pairs (Fig. 3) makes 
evident the reasoning of Hencken, Kuraev and SerboQ 
that unitarity corrections are small for e + e _ cross sec- 
tions and large for Also clearly the region of nu- 
clear collisions (dashed line in Fig. 1) would provide no 
reduction of the e + e~ pair cross section (Fig. 2), but 
would slightly reduce the \i pair cross section (Fig. 1). It 
is obvious that even with a momentum dependent form 
factor there is significant contribution to the fi pair cross 
section here in the region of ion-ion overlap. This contri- 
bution must be eliminated for ultraperipheral collisions, 
and leads into the discussion in the rest of this section. 

As noted above, what states are considered in the uni- 
tarity consideration can be determined or defined by the 



energy scale of the detected particles. For example, when 
considering pair production, one might not con- 

sider the dominant soft e + e~ pairs as part of the the 
excited spectrum for purposes of unitarity normaliza- 
tion. This might be a reasonable definition corresponding 
to the experimental detection conditions. On the other 
hand, to construct a calculated cross section correspond- 
ing to the observed pair production events without any 
other final state particles arising from ion-ion overlap, 
one should include only ultraperipheral impact parame- 
ters. 

In terms of a b (impact parameter) dependent ampli- 
tude M(k, b) exp(zk • b) an appropriate non-unitarizcd 
probability can be written 



Pi{b) 



d 2 fcM,(k,b)exp(ik-b)| 2 



(14) 



Let Pi (6) be a non-unitarized probability for exciting a 
pair and Pj(b) the corresponding probability for a 
nuclear reaction. Then define a partially exclusive cross 
section as one that excludes nuclear interaction only 



d 2 bP?(b) 



d 2 bPi(b)exp(-P 3 (b)). (15) 



Since the nuclear interactions occur only below some 60 = 
i?i + i?2, is convenient to express this cross section as a 
difference 



d 2 bP t (b) 



b<b 



d 2 bPi(b)(exp(Pj(b) - 1)). 



(16) 



The first term can be evaluated as was done before with- 
out a specific impact parameter representation 18]. The 
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TABLE I: RHIC: Au + Au, 7 = 100, fi + fi" total cross section. 



RESULTS IN mb 


Perturb. Exact 


b independent formulation 
b integration (ion overlap) 
b independent minus ion overlap 
Hencken et al. (7 = 108) 


211 164 
232 (36) 181 (20) 
175 144 
230 230 



TABLE II: LHC: Pb + Pb, 7 = 2760, ^+ \T total cross sec- 
tion. 



RESULTS IN barns 


Perturb. Exact 


b independent formulation 
b integration ion overlap 
b independent minus ion overlap 
Hencken et al. (7 = 3000) 


2.43 2.09 
0.13 0.06 
2.29 2.03 
2.60 2.60 



second term can then be evaluated using the method of 
Ref. [l!|: since b is limited to the lowest impact param- 
eters the exp(ik • b) factor is still numerically tractable 
even though k scales up by a factor of m^/m e as com- 
pared to the e + e~ case. A similar trick to subtract the 
nuclear interaction at small impact parameters has pre- 
viously been used [13] • 



V. NUMERICAL RESULTS 



cr(e + e~) from the perturbation theory result fl8j . Due 
to the higher values of transverse momentum transferred 
from the virtual photons in this LHC case it was not fea- 
sible to compute the b-dependent cross section contri- 
butions throughout the entire impact parameter range. 
However in the region of ion overlap the impact param- 
eter was small enough that the the a rapidly oscillat- 
ing phase exp(ik • b) in the integral over the transverse 
momentum k transfered from the ion to the lepton pair 
remained tractable and Eq. (16) could be utilized. The 
additional reduction from exclusion of overlap impact pa- 
rameters was 0.06 b for a best value of 2.03 b, an overall 
16% reduction from perturbation theory. 

The present perturbative a(p + p~) calculations are in 
fair agreement with the calculations of Hencken, Ku- 
raev, and Serbo[4(, but the present exact cross section 
calculations are in disagreement with their argument 
that Coulomb corrections are relatively insignificant for 
fi pairs. In this work I have shown that unlike the case 
for e + e~ pair production, the finite size of the colliding 
nuclei provides an important modification for both the 
perturbative and higher order calculated total cross sec- 
tions a{n + fj,~). The form factor reduces the higher order 
calculation by an even greater percentage than it does for 
perturbation theory. Furthermore, making the necessary 
elimination of interactions where the ions overlap further 
reduces the higher order cross section from perturbation 
theory. 



Table I summarizes the total cross section results for 
pair production at RHIC. In the first row the b- 
independent perturbative cross section for a(fi + fi~) is 
211 mb. and the cross section with higher order effects 
is 164 mb. As previously noted in Section III, this re- 
duction of 22% from perturbation theory is greater than 
the 17% reduction in the exact tr(e + e~) from the per- 
turbation theory seen in Ref. 18]. The second row shows 
the results from integrating the b-dependent computa- 
tion of the cross section shown in Fig. 2. The numbers 
in parentheses correspond to the negative of the second 
right hand term of Eq. (16). A rough check can be 
done by comparing the results of the b-dependent and 
b-independent calculations of the total production 
cross section. Even though the b-dependent becomes 
more inaccurate beyond the low impact parameters it 
still reproduces the b-independent results to about 10%. 
The lowest impact parameters have the greatest accuracy 
and we calculate the higher order cross section from the 
overlap impact parameters to be subtracted off as 20 mb. 
Thus the final computed best cross section (third row) is 
144 mb, a 32% reduction from the perturbation theory 
calculation. 

Table II shows calculations of colliding Pb + Pb ions 
at the LHC. The perturbative production cross 

section shown in the first row is 2.43 b and higher or- 
der effects reduced it by 14% to 2.09 b. Again this re- 
duction is greater than the 11% reduction in the exact 
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APPENDIX A: HIGHER ORDER FORM FACTOR 
EFFECTS 

As noted in Section III, to include a form factor in 
the eikonalized expression for the transverse integral with 
Coulomb corrections, the most obvious prescription is to 
apply the form factor to the transverse potential, leading 
to Eq. (9): 

Jtf(k) = 2tt f dppJ (kp){exp[2iZag(k)Ka(pu;/ 7 )} - 1}. 

( A1 ) 

But including the form factor in the transverse poten- 
tial is equivalent to letting the coupling constant Za run 
as a function of k, analogous to the situation in QCD. 
To do this makes the numerical integration of Eq. (1) 
more complicated and has not been done in this paper. 
However, by a relatively simple modification of the b in- 
dependent expression for the higher order cross section 
one can put an upper limit on the modification to higher 
order effect of using the more proper Eq. (9) rather than 
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the expression Eq. (10) utilized throughout this this pa- 
per. 

The numerical integration of Eq. (1) is most conve- 
niently carried out after a change of variables to £ = kp 
and 

F(k) = |J J ^J ^){exp[2iZaK (^/ 7 k)] - 1}. 

This integral is carried out for various values of the pa- 
rameter kj/ui. F actually has a two dimensional pa- 
rameterization in kj/u and k, but the 1/k 2 dependence 
trivially factors out of the integral. Likewise the addi- 
tional k dependence in expression Eq. (10) utilized in 
this paper factors out trivially. However the additional 
k dependence in the more exact expression Eq. (9) does 
not factor out trivially, leading to the additional compli- 
cation. 

To put a limit on the error incurred by using the ex- 
pression Eq. (9) rather than Eq. (10), I begin by recall- 
ing that the organization of the ^-independent computer 
code utilized in Ref. [l8| involves a difference 

|AF(k) 2 | = |f(k)| 2 -|F (k)| 2 (A3) 

between the squared value of the higher order expression 
Eq. (1) and that of the perturbative expression, Eq. (2). 
Since the effect of the form factor g(k) is the same as 
reducing the value of Z as a function of k at large k (like 
running coupling in QCD) then evaluation of Eq. (A2) 
for various values of Z may be used as a proxy for the 
higher order dependence on g{k). That is, for a given Z if 
the form factor is reduced from unity by some percentage 
then it is equivalent to no form factor and just reducing 
Z by the same percentage. Numerical calculations of 
|AF^k) 2 | show that it scales as Z 4 for low values of Z 
and a little less than Z 4 as Z is increased. This scaling 
is consistent with the integral over (A3) 

G = J -^k 2 [\F{k)\ 2 -\F^)\ 2 ] (A4) 

in Lee and Milstein's analysis of higher order Coulomb 
corrections [lM [l6j . which takes the analytical form 

G = -%n{Za) 2 [Re^j{{l + iZa) + j B uler], (A5) 

where ^((1 + iZa) is the dig amma function and ^Euier 
is Eulcr's constant. This expression may be alternatively 
expressed as 

G = -8Tr{Za) 2 f{Za), (A6) 

where f(Za) is the same function that was presented by 
Bethe, Maximon and Davies[3] for Coulomb corrections 
to e + e~ photoproduction on heavy nuclei and takes the 
form 

oo 1 

f(Za) = {Zaf 7 2 + (7 Wv (A7) 
^ n(n z + (ZaV) 

G obviously scales as Z for low values of Z and a little 
less than Z 4 as Z is increased. 



The expression Eq. (A3) modified with the lowest or- 
der implementation of the form factor utilized in this 
paper takes the form 

|AF/°(k) 2 | = |F/°(k)| 2 -|F /0 (k)| 2 

= . 9 (fc) 2 (|nk)| 2 -|F (k)| 2 ). (A8) 

If one consider the analogous expression with the form 
factor to higher order, then replacing the g{k) 2 depen- 
dence of Eq. (A7) with g(k) 4 suggested by the Z 4 scaling 
seen in the difference without a form factor, 

|A^(k) 2 | = |i^(k)| 2 -|itf(k)| 2 

= g(k) 4 (\F(k)\ 2 -\F (k)\ 2 ), (A9) 

should slightly overstate the higher order effect of the 
form factor in Coulomb corrections. 

Recalculation of the exact b independent RHIC Au + 
Au cross section of Table I with the g{k) 4 form factor 
scaling of Eq. (A9) gives 171 mb, a 19% reduction from 
perturbation theory in comparison with the 164 mb 22% 
reduction using the more approximate g(k) 2 of Eq. (A8). 
Likewise for Pb + Pb at LHC the exact calculation with 
g(k) 4 gives 2.12 barns, a 13% reduction from perturba- 
tion theory in comparison with the 2.09 barn 14% reduc- 
tion with g{k) 2 . 

Both recalculations make only a small change from the 
lowest order treatment of the form factor in this paper. 
And since it is far from trivial to implement the more 
proper higher order treatment of the form factor of Eq. 
(9), especially in b dependent calculations, I have not 
done so in this paper. It seems that once the k depen- 
dent cutoff of the form factor is put in, then sharpen- 
ing the cutoff by an additional squaring has a relatively 
small effect. Even with the higher order g(k) 4 scaling 
calculations, the reduction from perturbation theory in 
a(p + p~) are still larger than the reductions in the anal- 
ogous exact er(e + e~) without a form factor from pertur- 
bation theory. 



APPENDIX B: LONGITUDINAL FORM FACTOR 
EFFECTS 

One might include longitudinal form factor effects by 
modifying Eq. (3) to make g(k) a function of k 2 + uj 2 /-f 2 
rather than simply a function of k 2 : 

9(k) = i+^+Lv^/A 2 - (Bi) 

I have recalculated b independent cross sections using 
Eqns. (5), (10) and (Bl) in place of (3), and I find a 
5% reduction for both the perturbative and higher or- 
der computations for RHIC but only a corresponding 1% 
reduction for LHC. The 5% reduction for RHIC is equiv- 
alent to calculations without a longitudinal factor, but 
with the value of A reduced from 80 MeV to 75.5 MeV. 
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